CESARO SUMMABILITY OF FOURIER ORTHOGONAL 
EXPANSIONS ON THE CYLINDER 
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Abstract. A result concerning the Cesaro summability of the Fourier or- 
thogonal expansion of a function on the cylinder, where the orthogonal basis 
consists of orthogonal polynomials, in the U' norms is presented. An upper 
bound for critical index S is obtained. 



1. Introduction 

Cesaro summability of an infinite series is a classical topic which still receives a 
considerable amount of attention. In particular, research into the Cesaro summabil- 
ity of Fourier expansions in a multivariate settingis quite active. For example, the 
Cesaro summability of a Fourier orthogonal expansion in orthogonal polynomials 
on the parabolic biangle is studied in [I], and the unit ball and sphere are studied 
in [2]. Investigations into approximation on the cylinder have also received some 
attention recently. In [3], approximation using Fekete and Leja points in studied, 
while in [9], an approximation technique relating orthogonal polynomials and the 
Radon transform on parallel disks was investigated. In this paper, we present a 
result on the convergence of the Cesaro means of a Fourier orthogonal expansion 
in orthogonal polynomials of a function defined on the cylinder. This result relies 
on results proven by Xu and Li in f6^, which concerns the Cesaro summability of 
functions on the hypercube. 

2. Background 

We denote the d-dimensional unit ball by B'^ and the m-dimensional hypercube 
[—1,1]™ by 7™. We first present the orthogonal polynomial basis used to obtain 
the Fourier orthogonal decomposition on the cylinder. We will use the notation 
Ud to denote the dimension of the space of orthogonal polynomials of degree n on 
either B'^ or which is known to satisfy 

nd = 

V n 

On the hypercube /"', the product Jacobi polynomials are used to form a basis. 
The univariate Jacobi polynomials Pn°"^\x) form an orthogonal polynomial basis 
on [—1, 1] with respect to the weight function w^^'^^x) = (1 — a;)"(l + x)^ , with 
a, (3 > —1; that is, 

(2.1) Pjr'^\x)Pj^^P\x)w^^'^\x)dx = 6^,rrCn,c,,p, 
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where 



2°+'3+i r(?i + a + l)r(?i + /3 + 1) 
2n + a + (3+lT{n+ l)r(n + a + /3 + 1) ' 



see [S] for more details. The orthonormal Jacobi polynomials are obtained by 
scaling the Jacobi polynomials so that the right side of (|2.ip is 6m,n', we denote 
the resulting orthonormal polynomials by pi^'^\x). Next, an orthonormal product 
Jacobi polynomial on /™ is defined by 



(2.2) 



3("-^^)(2;i)p^r'^^Ha;2) 



where we let a = (ai, a2, • • . , cUm), and similarly define /3, 7, and x. We will consider 
the degree of Pj°''^\x) to be the total degree, I7I := 71 + 72 + • ■ ■ + 7m- With the 
weight function w*^"'^)(x) — Hi" 1 '"^^"''^'■'(^^i); orthonormal basis for the space 
of polynomials of degree n on 7™ is obtained by taking the set of all polynomials 
of the form in (|2.2p . with I7I = n, where orthonormality is in the sense that 



(2.3) 



P^"'^)(x)P,("'^)(x)w("'^)(x)dx = S^^r,- 



It will be convenient to follow the notation developed in Chapter 2 of [1] and list 
the orthonormal basis of degree n polynomials on 7™ in column vector form. We 
define 



d(",/9) 
d(",/3) 



(x) 



(x) 



where 7"'^, 7^, . . . , 7"™ are multi-indices with I7I = n. 

On the d-dimensional unit ball B'^, we consider the weight function 'w^{y) = 
(1 — ||y|P)'^~^^^, with fi > 0, and denote a basis of orthonormal polynomials of 
degree n = \a\ on B"^ by s^i(y), s^a (y), ■ ■ • , s^^n^ (y), where orthonormality in a 
similar manner as (|2.3p . Several examples of specific orthonormal bases exist, see 
[Hp. 38]. If d = 1 and /j, = 0, these polynomials correspond to the Chebyshev 
polynomials of the first kind, while if d = 2 and /i = 1/2, an orthonormal basis is 
given by the polynomials 

2„/2-, + lp(0,«-2,)(2||y|p-l)5,,„_2^.(y), 

where < j < n and S';3,n-2j"(y) is a spherical harmonic on S"^, defined by 



:sin(n(7r/2-0)), 52,„(0) 



:Cos(n(7r/2-6')), 



with X = cos 9 and y = sin0; see [7] for more information on spherical harmonics. 
One important property of orthogonal polynomials on B'^ with respect to is the 
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compact formula introduced in [TO] . 



n + ^i+^ (r(/i)) 



(2.4) my)VW) = 



±_i 22M-ir(2/i) 



for > 0, and 



n 4- ^-J- r / d-i^ 



(2.5) [§^(y)] §°(y') = ^^r^ [ci^'(x • y + y^T^^T^) 

2 

(x • y - yr^b^yr^Bf ) 

for /i = 0. 

Following the column vector notation used on J™, we will use the notation 



(y) 



to denote the column vector whose elements form an orthonormal basis of polyno- 
mials of degree n on B'^ with respect to Wfj.. 

The reproducing kernel of degree n with respect to a weight function w on some 
set X of positive Borel measure, ii'„(x, x'), is a function satisfying 



/ /(x)/^„(x,x')«i(x)dx = /(x') 



when / is a polynomial of degree less than or equal to n. Given a basis of orthonor- 
mal polynomials with respect to w, {Pa{^)}aeA, for the space of polynomials of 
degree less than or equal to n, the reproducing kernel has the form 



/C„(x,x') = ^Fc(x)P„(x') 



The n'th partial sum of the Fourier orthogonal expansion of an integrable function 
/ on X with respect to w, S{w; /), is defined by 

5„(w;/)(x')= / /(x)i^„(x,x')?«(x)dx. 
Jx 

With the column vector notation above, the reproducing kernel on the hypercube 
with respect to w^"''^^ can be written as 



n 

if„(x,x')=^[pi-^)(x) 



k=0 



and a similar formula holds for the reproducing kernel on the ball with respect to 
Wfj,. We may then write the reproducing kernel on B'^ x /™ with respect to the 
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weight function w{a, (3, n;x,y) := w^"''^) (x)w^(y), i4:„(x, x', y, y'), as 

71 k 

i^„(x,x',y,y') = EEh"''^W 



'^Pf-'^^lx') k ,(y)l^§t,(y'). 



k=o j=0 

We define the n'th partial sum of the Fourier orthogonal expansion of an integrable 
function f on x /"* to be 

S'„(/U,Q;,/3;/)(x',y') = / / Jf^Cx, x', y, y')/(x, y)M;(a, ^, /z; x, y) rfxrfy. 

Jgd Jim 

Given a series X)s„, the Cesaro means, or (C,d) means, of the series is defined 
to be 

n j 
j=0 fe=0 

where cf^ j = -^z^^^, and {n)j = ni=i(^ + A; — 1). If we define 



Pf'«(x') St,(y) §^_,(y'), 



if;^(x,x',y,y')=X:<,5:[p^«(x)" 

fe=0 j=0 

then the Cesaro means of of order 6, or the (C, d) means, of the Fourier orthogonal 
expansion of / are defined by 

5^(M,a,/3;/)(x',y')= / / <(x,x',y,y')/(x,yHa,^,/i;x,y)dxrfy. 

We will be investigating the value of 6 for which this series converges in the space 
LP{B'^ X I"^;w{a, I3fj,;x,y)) - that is, for what 6 is 

lim / |S'^(a,^,/i;/)(x,y) - /(x,y)ru;(a,^,/i;x,y)dxrfy = 0. 

3. Main Theorem 

Wc now present our main result for this paper. 

Theorem 3.1. Let f be a continuous function on B'^ x I™, and suppose that 
H > 0, ai > —1, j3i > —1, and ai + Pi > —1 for 1 < i < m. The Cesaro means 

of the Fourier orthogonal expansion of f with respect to w{a, /3, ii;x,y) converge in 
LP{B'^ X I"';w{a,P,ii;yi,y)), with 1 <p < oo, and C{B'^ x J™), to f if 



(3.1) 5 > ^max{ai,/3i} + 



d + m—1 



i=l 



+ max s 0, — ^ min {aj, — jj, ■ 



d + m+1 



Proof. Our proof will ultimately reduce the question of convergence on the cylinder 
to that of the hypercube. We will show that the integral 

/ / |-^n(x,x',y,y')| w(a,/3,/x;x,y)dxrfy 

J/™ JSd 

is uniformly bounded by some constant M which is independent of n, x', and y', 
which will then imply summability, by results known on the hypercube. Throughout 
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our proof, c will denote a positive constant that may change values from line to 
line. 

Following Lemma 2.2 in we first show that it is enough to consider x' = e := 
(1, 1, . . . , 1) in the hypercube. 

Lemma 3.2. In order to prove the convergence of {C,d) means of the orthogonal 
expansion, it suffices to prove 

(3.2) / / |A'^(x,e,y,y')|«i(a,/3,/x;x,y)dydx<M 

for M independent of n and y' . 

Proof. The proof of the lemma follows from the following theorem, which appears 
in [51 p. 262]. 

Theorem 3.3. Let a,/? > —1 and a> (3. An integral representation of the form 
PlC'P\x)PlC^P\y) = C P,("'«(l)P^"''3)(z)i^("''3)(a;,y,z)u;("'^)(^)'^^ 



exists, with the function K satisfying 

K'^^^f^^ (x, y, z) w^"-'^^ (z) dz < M 

J -I 

for —1 < x,y < 1. 

This result is for univariate Jacobi polynomials, but easily extends to the product 
Jacobi polynomials as 



Pi"'« (x)Pi"'^) (y) = j^^^ P^"'« (e)Pi"'« (z)K("'« (x, y, z)«;("'« (z)dz, 
where x, y G F^, and 

in 

K("'«(x,y,z)=ni^^"-''''^(^^, 2/^,^0, 
1=1 

and K("'«(-,-,-) satisfies 

|K("^'3)(x,y,z)|w("^'3)(z)dz < M, 

where M is a constant given by the product of the constants in the univariate case. 
We obtain 



^« — / / |^fi(x,x',y,y')| w(a,^,Ai;x,y)dydx 



j=0 fc=0 

X w{a, /3, ^; , x, y) dydx 



< 



/™ JB<i J I 



Pj^f^Cz) [§^(y)]^§t'(y') 



j=0 k=0 

X |K(x, x', z)| w''°''^\z)dz w{a, jS, ^jl; x, y) dy dx. 
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Applying Fubini's theorem gives 

n j 



At < 



j=0 fc=0 



X J |K(x, x', z)| w^"'^^(x) dx w(a, /3, /i; z, y) dz dy 
<M [ [ \Kf^{e,z,y,y')\wia,/3,ii;z,y)dydz. 



Replacing z with x and switching the places of e and x proves the lemma. 



□ 



Our next lemma reduces the integral over B"^ to an integral over [—1,1] of a 
Gegenbauer polynomial. The idea for this lemma comes from the proof of Theorem 
5.3 in p]. We define 

Qlr^^\y')--^ I K(x,e,y,y')|z«^(y)dy 



and 



j=0 fe=0 ^ + 2 

Lemma 3.4. For /i > 0, 

»i 



(3.3) 



Proof. We first consider the case /i > 0. Substitute (|2.4p into (|3.2I) to obtain 



(3.4) (?("■«(/) 



Fi ((y, y') + v/l-|y|Vl-|y'^) (1 - dt 



Applying the change of variable y = rrj, where rj E S"^ ^,0<r<l, gives 

r-l 



Sd-i 



\ Ft {r{ri.y') + ^l-\yf^/l^t^ (1 -i^Y ' dt 



(l-r2)'^-i/2dcj(7y) dr, 



where doj is the surface measure on S''^ ^. Now let A be the rotation matrix satis- 
fying A{y') — (0, 0, . . . , 0, |y'|), and apply the change of basis rj i-7> A^-q to obtain 



e^''^y')- / 



I r'-' J^^_^ J ^ Fi [rr^a W\ + ^1 - |y|' v/T^i) (l - 



dt 



X (l-r2)^-i/2^cj(?7)dr 
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where 77 = (771, ... , r]d)- If we let r]d — s, then 77 = (Vl — 5^7, s) for some 7 e S''^ ^, 
and changing variables gives 

(3.5) ^^''Hy')=^d-2_^'r'-'/_'^ 

X (^rs \y'\ + y^l - |yf ^137^^^ (1 _ ^2)^-1 

where oJd-2 is the surface area of S'''"^. Let s n- p/r so ds = dp/r and move the 
absolute value inside the innermost integral to obtain 



Switching the order of integration of r and p and applying the change of variable 
q Vl — r'^t, dq = \/l — rft gives 



1 — r —q) dqr[r —p)^~drdp. 
Switching the order of integration of q and r gives 



li-Vi-W 



i^„^(^p|y'| + Vl-|y'l'9 

/ /1 2 2\'*~1 / 2 2 

/ (1 — r — (7 j r(r — p 

J\p\ 



2 dr dp. 



Applying the change of variable = u(l — q^ — p^) + p^ shows the inner integral 

is 5(1 — q^ —p^Y'^ 2 B{ii, ^^), were -B(x, y) is the beta function. Hence, we have 
the inequality 



Wd-2 



iJ-y/i-\p? 



Fnpw\ + ^i-\yrq 



X (l-g2_p2^- 



dqdp. 



Next, we apply the change of variable g — p^s to obtain 

(3.6) t/("'^)(y') < c J'^ (p |y'| + ^l-|yf yr^s) 

d — 2 ■ d — 3 I 

X (l-p2)^+'^ (l_s2)^+'' 



dsdp. 
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Changing variables once again, we let u = p |y'| + y 1 — |y'|^^l — p^s to obtain 
(3.7) e^^^y')<c/ / H{u)\D^.{\y'lp,u) 



1 "'piy'i-Vi-iy'lVi-p^ 



2\ — — I-/' 



du dp, 



where the function D\(v,p,u), introduced in 10 , is defined by 



A-l/2 



Dx (t',p, u) = 



(l — i;^ — — + 2upv) 



[{l-v^)il-u^){l-p^)Y 

for 1 — ti^ — — + 2upti > and otherwise. It is readily verified that 
»i 

Dx {u, v,p) (1 - dp ^ 2^^B{X + 1/2, A + 1/2). 
Hence switching the order of integration in p.7p . we have 
(3.8) gl"'^Hy')<c [' \FfAu)\ f D^^(\y'\,u,p){l-p')'^^' dp 



-1 

X fl — u 
1 



d-2 I 

2\— — l-Ai 



du 



<c \FfM\{l-u') 



du. 



This proves the lemma for /i > 0. 

Turning our attention now to the case when /i = 0, we substitute (|2.5p into the 
left side of p.2p and ignore the integral over I™ as before to obtain 



gr\y') = I 

JB' 



i-|ylVi-ly'l' 



i^,? (y,y')Vi-|yrVi-|y' 



wo{y)dy. 



We perform the same change of variables from the case when > to obtain the 
equivalent of p.Sp . 







2 f'r'-' 








Jo .J 





F' 



(rs\y'\ + Vl-|y'|V 



1 — r' 



+ F' 



{l-rY^^^(.^-s')~ dsdr. 



d-3 



Now we substitute p = \/l — and let w = ^/l — |y'p to obtain 



„ (i-2 „ <i-3 

(1 -p2)~(l - s2)~ dpds 



(3.9) 



„ d-2 „ <i-3 

X (1 -p^) 2 (1 - s^) 2 dpds. 
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The right side of p.9|) is the right side of p.6p . with v in place of |y'|. FoUowing 
the same steps of the proof for /i > 0, we obtain the equivalent of (|3.8|) . 



□ 



g[r'^\v')<cj\Fi{u)\{\-u')^du, 

which proves the case for = 0. 

To finish the proof of Theorem 13. 1[ we substitute p.3p into 
/ / |^n(x,e,y,y')| w(a,/3,/i;x,y)dy dx 



to obtain 



<c 



"^1" X! 

j=0 fc=0 



A* ^ 2 



[Pj-f) (e) 



[l-u 



d-2 
2N — 



<iMw("^^)(x) dx. 



After substituting in the well known identity for Gegenbauer polynomials, 



n + \ 
A 



we arrive at the following inequality. 

/ / |-^*(x,e,y,y')| w7(a,/3,/i;x,y) 



<c 



j=0 k=0 



H-f^ (x) 



Fj-f ) (e) 



(1-^2)—+^ d«w(°^^)(x)d> 



Since the Gegenbauer polynomials are a subset of the Jacobi polynomials, the proof 
of our theorem then follows from Theorem 1.1 and Lemma 2.2 from [5], which , 
when combined, can be restated as follows. 

Theorem 3.5. Let a,; > —1, /3i > —1, and at + Pi > —1 for 1 < i < m. If 

m + 2 1 



5 > y^max{ai,/3J 

i=l 

then the integral 



],(",/3) 



(x) 



U'("''3)(x)rfx. 



is uniformly bounded for all n and x G which is sufficient to show that the 
Cesdro means of the Fourier orthogonal expansion of a continuous function f on 
with respect to w^"'*^) converge in the norms of LP{I"^;w^°''^'^ ) for 1 < p < oo 
and Cil""). 



□ 
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